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ABSTRACT 

We propose the Galactic flare model for the origin of the X-ray gas in the Galactic halo. For 
this purpose, we examine the magnetic reconnection triggered by Parker instability (magnetic 
buoyancy instability), by performing the two-dimensional resistive numerical magnetohydrody- 
namic simulations. As a result of numerical simulations, the system evolves as following phases: 
Parker instability occurs in the Galactic disk. In the nonlinear phase of Parker instability, 
the magnetic loop inflates from the Galactic disk into the Galactic halo, and collides with the 
anti-parallel magnetic field, so that the current sheets are created in the Galactic halo. The 
tearing instability occurs, and creates the plasmoids (magnetic islands). Just after the plasmoid 
ejection, further current-sheet thinning occurs in the sheet, and the anomalous resistivity sets 
in. Petschek reconnection starts, and heats the gas quickly in the Galactic halo. It also creates 
the slow and fast shock regions in the Galactic halo. The magnetic field (B ~ 3 fxG), for 
example, can heat the gas (n ~ 10~ 3 cm -3 ) to temperature of ~ 10 6 K via the reconnection in 
the Galactic halo. The gas is accelerated to Alfven velocity (~ 300 km s _1 ). Such high velocity 
jets are the evidence of the Galactic flare model we present in this paper, if the Doppler shift 
of the bipolar jet is detected in the Galactic halo. 



Subject headings: Galaxy: halo — ISM: magnetic fields 
magnetohydrodynamics 



instabilities 



1. INTRODUCTION 

The X-rays from hot g observed in the Galac- 
tic halo. Its luminosity and temperature are Lx ~ 
7 x 10 39 erg s" 1 and T ~ 10 6 K (Pietz et al. 1998). 
The volume and thermal energy are estimated to be 
E ~ 10 55 erg and V ~ 10 68 cm 3 . To explain such 
hot gas in the Galactic halo, the "Galactic foun- 
tains" model has been proposed (i.e., supernova rem- 
nants and stellar winds heat the gas; Bregman 1980; 
Norman & Ikeuchi 1989; Shapiro & Benjamin 1991; 
Shull 1996; de Avillez 2000; Slavin, McKee, k Hol- 
lenbach 2000). The energy source, however, may not 



be explained fully by this model because the evidence 
is not observed adequately (see also Birk, Lesch, & 
Neukirch 1998). In this paper, we propose another 
mechanism. It is the Galactic flare model, i.e., the 
magnetic heating in the Galactic halo. 

Parker (1992) pointed out the importance of mag- 
netic reconnection for the heating of Galactic plas- 
mas. Whenever the magnetic flux collides with an- 
other flux, which is not exactly parallel, the current 
density increases, and a strong dissipation sets in 
(e.g., via anomalous resistivity) to trigger fast recon- 
nection (e.g., Ugai 1986). The magnetic reconnec- 
tion is observed in solar flares by the X-ray satellites 
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Yohkoh (Masuda et al. 1994; Shibata 1996; Tsuneta 
1996) and SoHO (Yokoyama et al. 2001). In the solar 
flare, the reconnection heats the plasma from a tem- 
perature of ~ several xlO 6 K to ~ several xl0 7 ~ 8 
K, and accelerates it to Alfven velocity (~ 10 2 ~ 3 
km s _1 ) (e.g., Shibata 1996). The reconnection 
would occur also in the Galaxy (Tanuma et al. 1999a, 
1999b, 2001a, 1999b), which may be called "Galactic 
flare" (Sturrock k Stern 1980; Kahn k Brett 1993). 
Total magnetic energy is E mag ~ ((5) 2 bs /87r) V G 
~ 10 54 ' 4 erg at least, where (.B)obs(~ 3 fiG) is the 
mean observed field strength (see Beck et al. 1996; 
Vallee 1997), and V G (~ 10 67 cm 3 ) is the volume 
of the Galaxy. The rotational energy of the Galaxy 
(~ io 58 - 9 erg) and kinetic energy of the interstellar 
gas (~ 

10 58.2 

erg) are its origin (e.g., Parker 1971; 
Sturrock k Stern 1980; Tanuma et al. 1999a; Tanuma 
2000). 

The steady reconnection mechanisms were pro- 
posed (see Priest k Forbes 2000). In Sweet(1958)- 
Parker(1957) type reconnection, the diffusion region 
is so long as to occupy whole current system. It 
can not be applied to the solar flare phenomena, be- 
cause the reconnection rate of this model is too small 

— 1/2 

(~ R m ) in the solar corona, where R m is the mag- 
netic Reynolds number (~ 10 12 ). On the other hand, 
in Petschek(1964) type reconnection, the diffusion re- 
gion is localized near an X-point, and standing slow 
shocks occupy whole current systems. In this case, 
the energy conversion via slow shocks is much larger 
than Ohmic heating. This can hence be applicable to 
the solar flare phenomena, because the reconnection 
rate of this model is ~ 0.1 — 0.01. This is called 
"fast reconnection". A basic problem of Petschek 
model is as follows: Petschek reconnection occurs, 
if the anomalous resistivity sets in the current sheet 
(e.g., Ugai 1986; Yokoyama k Shibata 1994; Tanuma 
2000; Tanuma et al. 1999a, 2001a). The anomalous 
resistivity set in, when the current-sheet thickness be- 
comes comparable with ion Lamor radius or ion in- 
ertial radius. It is, however, not fully known how 
the current sheet becomes thin, because the typical 
size of solar flare (10 9-11 cm) is much larger than 
these radii (10 2-3 cm). This situation is similar to 
that of the Galaxy (R m > 10 15 ), where typical size of 
magnetic field (> 10 19 cm) is much larger than the 
ion Lamor radius (~ 10 7 cm). To solve these prob- 
lems, we proposed the current-sheet thinning via the 
"fractal tearing instability" (Tanuma 2000; Shibata 
k Tanuma 2001; Tanuma et al. 2001a). 

Many two-dimensional (2D) magnetohydrody- 
namic (MHD) numerical simulations have been car- 
ried out for the magnetic reconnection in the solar 
atmosphere (Magara k Shibata 1997; Odstrcil k Kar- 
licky 1997), and in the Galactic halo (Zimmer, Lesch, 



k Birk 1997; Birk et al. 1998), by assuming the cur- 
rent sheet at the initial condition (see also Nitta et al. 
2001). Recently, Tanuma et al. (1999a, 1999b, 2001a) 
examined the magnetic reconnection triggered by a 
supernova-shock by performing the 2D MHD simu- 
lations with a high spatial resolution, and proposed 
that it can generate X-ray gas in the Galactic disk 
(e.g., Ebisawa et al. 2001). They found that the tear- 
ing instability (Furth, Killeen, k Rosenbluth 1963) 
occurs in the current sheet long after the passage 
of a shock wave. Petschek reconnection occurs after 
further current sheet thinning via secondary tearing 
instability. In the present model, Parker(1966) in- 
stability creates the current sheet by itself and trig- 
ger the magnetic reconnection (Tanuma et al. 1999b; 
Tanuma 2000; see also Shibata, Nozawa, k Mat- 
sumoto 1992; Yokoyama k Shibata 1996, 1997). 

Recently, we examied three-dimensional (3D) MHD 
simulations of the magnetic reconnection with a low 
spatial resolution. Petschek reconnection occurs after 
the current sheet thinning by the tearing instability 
in both 2D (Tanuma et al. 1999b) and 3D models 
(Tanuma 2000; Tanuma et al. 2001b), because we 
can not resolve the secondary tearing instability when 
we assume a rough grid (The similarities between 2D 
and 3D models are consistent with Ugai k Shimizu 
1996). 3D effect such as Rayleigh- Taylor instability, 
however, appears when reconnection jet collides with 
high pressure gas and magnetic loop much after the 
onset of Petschek reconnection. Tanuma et al. (2002) 
applied the results to the creation of helical magnetic 
field and confinement of high energy particles in the 
solar flare. We study the basic physics of magnetic 
reconnection which are common between 2D and 3D 
models, although 2D model examined in this paper 
is a toy model of limited in 2D dimension. In this 
paper, however, we examine the 2D model under a 
higher spatial resolution than the 3D model which 
we are able to do. 

Parker instability is the undular mode (k \\ B) of 
magnetic buoyancy instability (Parker 1966), which 
occurs if a gas layer in a gravitational field is sup- 
ported by the horizontal magnetic fields. Suppose 
that the magnetic field lines are disturbed and be- 
gin to undulate. The gas in the loop top slides down 
along the field lines, so that loop rises further, and 
the instability sets in. Parker instability is suggested 
to influence the motion of clouds, HII regions, and 
OB associations (Tosa k Sofue 1974), and the dis- 
tribution of clouds (Mouschovias, Shu, k Woodward 
1974; Vrba 1977; Blitz k Shu 1980; Elmegreen 1982); 
for example, Perseus hump (Sofue k Tosa 1974), 
Perseus arm (Appenzeller 1971, 1974), Barnard loop 
(Mouschovias, Shu, k Woodward 1974), and Sofue- 
Handa(1984) lobe. 
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The linear analysis of Parker instability were made 
by many researchers (Shu 1974; Horiuchi et al. 1988; 
Hanawa et al. 1992; Foglizzo & Taggar 1994; Chou 
et al. 1997; Kamaya et al. 1997). The 2D MHD 
simulations were performed for solar flares (Kaisig 
et al. 1990; Nozawa et al. 1992; Shibata et al. 
1992; Yokoyama &, Shibata 1996), and the interstel- 
lar medium (Basu et al. 1997; Matsumoto et al. 1998; 
Santillan et al. 2000; Steinacker & Shchekinov 2001). 
The three-dimensional (3D) simulations of Parker in- 
stability of horizontal magnetic field in solar atmo- 
sphere and Galaxy (Matsumoto et al. 1993; Kim, 
Ryu, & Jones 2001; Hanasz, Otmianowska-Mazur, &: 
Lesch 2002), and the twisted flux tube in the solar 
atmosphere (Matsumoto et al. 1998; Abbett, Fisher, 
& Fan 2000; Fan 2001; Magara 2001), the Galaxy 
(Hanasz & Lesch 2000; France et al. 2002), and accre- 
tion disks (Ziegler 2001) are also performed. Shibata 
et al. (1989, 1992) and Yokoyama & Shibata (1994, 
1996, 1997) examined the magnetic reconnection trig- 
gered by Parker instability in the solar corona. Re- 
cently, Hanasz et al. (2002) examined the 3D model 
of the magnetic reconnection in magnetic loop cre- 
ated by Parker instability with Coriolis force in the 
interstellar medium. In the present paper, we extend 
Shibata et al. (1989, 1992) and Yokoyama & Shibata 
(1994, 1996, 1997) 's solar flare model to the Galaxy: 
The Galactic flare as the origin of the X-ray gas in 
the Galactic halo. 

In this paper, we propose a possible origin of X- 
ray gas in the Galactic halo. In the next section, we 
describe the simulation method. In sections 3 and 
4, we describe the results of numerical simulations, 
and discuss it. In the last section, we summarize this 
paper. 

2. NUMERICAL SIMULATIONS 

2.1. The Situation of the Problem 

Figure 1 shows our schematic scenario. Figure la 
displays the initial condition. Parker instability oc- 
curs in the Galactic disk (Fig. lb; we call this sit- 
uation phase I and II in this paper). The magnetic 
reconnection occurs and heats the gas (Fig. lc; phase 
III) . The heated gas is confined by the magnetic field 
(Fig. Id). 

2.2. Two-Dimensional Resistive MHD Basic 
Equations 

The resistive MHD basic equations are written as 
follows: 



dp 
dt 



+ V-(pv) = 0, 



(1) 



p-^ + p(v-V)v + V Pg = -JxB + pg, (2) 



dB 

— V x (v x B) 

dt y ' 

de _ 

_ + V-[(e + p>] 



-cV x ( V J), (3) 
r]\J\ 2 + v-V Pg (4) 



where p, v, B, n, e, J, and g are mass density, veloc- 
ity, magnetic field, electric resistivity, internal energy, 
current density (= cV x B/Air), and gravity, respec- 
tively. We use the equation of state for the ideal gas, 
i.e., p g = (7 — l)e, where 7 is the specific heat ratio 
(=5/3). We take Cartesian coordinates (x,z), where 
x and z axes are in horizontal and vertical directions, 
respectively. The velocity, magnetic field, and gravi- 
tational acceleration are v = (v x ,v z ), B = (B X ,B Z ), 
and g = [0, — <7tanh(z/0.001)], where g is constant. 

2.3. The Normalization 

We normalize the velocity, length, and time by the 
sound velocity (C s ), gravitational scale height (H), 
and H/C s , respectively. The units of normalizations 
are 



C s ~ 10T 4 1/2 km s~\ 
C 2 

H ~ s. ~ ioor! /2 P c, 
9 

H 7 
r = — ~ 10' yr, 



(5) 
(6) 

(7) 



where T4 is the temperature in unit of 10 4 K, and g is 
the acceleration of the gravity (Tanuma et al. 1999b; 
Tanuma 2000). The units of density, gas pressure, 
magnetic field strength, current density, and resistiv- 



ity are p 

erg cm" 3 , (p C 2 s /M) 1 / 2 



1 x 10 25 g cm 3 , p C, 



10 



-13 



3.2 pG, (poC^/S^/H, 
and c 2 HC s /4:TT, respectively. In this paper, we solve 
the non-dimensional basic equations, and describe the 
results of numerical simulations by non-dimensional 
values. 



2.4. The Anomalous Resistivity Model 

We assume the anomalous resistivity model as fol- 
lows: 



V 



if v d < v c 

a(v d /v c - l) 2 if v d > v c 



(8) 



(e.g., Ugai 1986, 1992; Shibata et al. 1992; Yokoyama 
& Shibata 1997; Ugai & Kondoh 2001; Tanuma et al. 
1999b; Tanuma 2000), where v d {= \J\/ p) and v c are 
the relative ion-electron drift velocity and threshold 
above which anomalous resistivity sets in. The pa- 
rameters are a = 0.01 and v c = 400. 
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2.5. The Initial Condition of Typical Model 

For the initial condition, we assume that the cool, 
dense Galactic disk is between the hot, rarefied 
Galactic halo (Fig. la; Horiuchi et al. 1988; Mat- 
sumoto et al. 1988; Tanuma et al. 1999b; Tanuma 
2000) . Tables 1 and 2 show the variables and param- 
eters. The temperature is 



T(x,z) = T + 0.5 



tanh 



z — z c 
0.5 



+ 1.0 



{T c ~ T 



(9) 

The parameters are To = 1, T c = 25, and z c = 6. 
The basic physics does not much depend on T c /Tq 
(Matsumoto et al. 1988; Yokoyama & Shibata 1996). 
The magnetic fields are in the Galactic disk {\z\ < 
Zd = 2.5) and in the Galactic halo (\z\ > z^ = 5). 
The initial magnetic field is horizontal (B z (x, z) = 0) 
everywhere, so that the ratio of gas to magnetic pres- 
sure is /?o = 8irp g (x,z)/B x (x,z) 2 = 0.2. The (3(z) 
is 



P(z) 



+ 



1 

Wo 
1 

Wo 



1 — tanh 



1 + tanh 



z 



Zd 



0.5 

z - z h 
0.5 



(10) 



The value {[5q = 0.2) is attained in the Galactic disk, 
if the magnetic field of B ~ 5 1 / 2 (S) bs ~ 7 pG, where 
(B) obs is the mean observed field strength (~ 3 pG) in 
the Galactic disk (see Beck et al. 1996; Vallee 1997). 
The magnetized gas is under the MHD equilibrium; 



d_ 

dz 



Pg{z) + 



B x (zf 
8vr 



+ p(z)g = 0. 



(11) 



The variables are po = 1, p g o = 0.6, and B x q ~ 8.68 
at the equatorial plane. The sound and Alfven ve- 
locity are C s = (jp g o/ Po) 1 ^ 2 = 1-0 and v 1 ^ = 
-Bxo/(4vrp ) 1/2 ^ 2.45 in the Galactic disk. The 
Galactic disk is unstable to Parker instability. 

We use the N x = 403 grids and N z = 604 grids 
in the horizontal and vertical directions, respectively. 
The intervals of grids are uniform (Ax = 0.30) in 
x-axis and nonuniform (Az > 0.075) in z-axis. We 
assume the top (z = +34.5) and bottom (z = —34.5) 
surfaces are free boundaries, and the right {x = 
+30.0) and left (x = —30.0) ones are periodic ones. 
We use the 2-steps modified Lax-Wendroff method. 
We put random perturbations (< 0.05C S ) on vertical 
velocities [v z (x,z)] in the Galactic disk. 

We neglect radiation and heat conduction. The 
cooling times due to the conduction and radiation 
for the cool (T ~ 10 4 K), dense (n ~ 0.1 cm~ 3 ) gas 



in the Galactic disk are 
nkTX 2 



''"cond 



K 
10 



T 7/2 



13 



n 



0.1 cm 3 
A \ 2 / T 



Trad 



100 pc 
nkT 
n 2 A(T) 



10 d 



10 4 K 



-5/2 



(12) 

vr. (13) 
(14) 



n 



10 4 KJ VO.l cm" 3 



A(10 4 K) 



10 -21 erg cm 3 s~ x 



yr, 



(15) 



respectively, where A(T) is the cooling function 
(Spitzer 1978), and kq is constant (= 10~ 6 erg s _1 
cm" 1 K" 1 ). For the hot (T ~ 10 6 K), rarefied 
(n ~ 10~ 3 cm -3 ) gas in the Galactic halo, however, 
they become 



7"cond 



10 a 



n 



10 3 cm 3 



A 



eff 



3kp C y V10 6 K 



-5/2 



(16) 



''"rad 



10 y 



n 



10 6 K 
A(10 6 K) 



10 -3 cm 3 
l 



10 -23 erg cm 3 s" 1 



yr, 



(17) 



respectively, where A e ff is the effective length of heli- 
cal magnetic loop. They are one order of magnitude 
longer than the typical time scale (~ 10 8 yr; see next 
section) of the physical process examined in this pa- 
per, so that the cooling mechanisms can be neglected, 
once they are heated to X-ray gas in the Galactic 
halo. The basic properties of the magnetic recon- 
nection such as reconnection rate and energy release 
rate are not much affected by the heat conduction 
(Yokoyama & Shibata 1997). 

3. THE RESULTS OF NUMERICAL SIMULATIONS 
3.1. The typical Model (Model Al) 

Phase I: The Linear Phase of Parker Insta- 
bility (t < 40) 

Figures 2-4 show the time evolution of the sys- 
tem. The axes are in the unit of ~ 10H ~ 1 
kpc. The magnetic field lines starts to bent across 
the plane. Figures 2a, 3a, and 4a shows a char- 
acteristic feature of the odd-mode (glide-reflection 
mode), which grows earlier than even- mode (mirror- 
symmetry mode) (Horiuchi et al. 1988; Matsumoto 
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et al. 1988; see also Tanuma et al. 1999b; Tanuma 
2000). The magnetic field inflates toward the Galac- 
tic halo by the magnetic buoyancy force. The gas 
slides down along the magnetic field lines. 

Phase II: The Nonlinear Phase of Parker In- 
stability (t ~ 40 - 60) 

The system enters the nonlinear phase at t ~ 40. 
In the valleys of the waving field, the vertical dense 
spurs are formed almost perpendicular to the Galac- 
tic plane (Fig. 3 a). The dense regions are also created 
in the valleys. Figure 5 shows the time variations of 
the drift velocity and velocity. The velocity increases 
to Alfven velocity (~ 2.5 — 3.0) in this phase. 

Figure 6 shows the time variation of the various en- 
ergies. The magnetic, thermal, kinetic, gravitational, 
and total energies are defined by 



where n m , n out , and T m are the densities of inflowing 
and outflowing gas, and the temperature of inflow- 
ing gas. It is derived from the equation 
n- m kT m + B 2 /8ir. We assume (3 ~ (3q ~ 0.2, n m ~ 



n out , and T u 



25. 



Figures 7 and 8 show the slow and fast shock re- 
gions, respectively. Two peaks in the profile of cur- 
rent density are a characteristic pattern of the slow 
shock (see also Fig. 4c). The fast shock is also cre- 
ated by the collision between the reconnection jet and 
ambient high-pressure gas (Fig. 8). 

The magnetic energy is released at the rate of 
—dE m&jg / dt ~ 35 (Fig. 6). It is determined by the 
Poynting flux entering into the reconnection region, 



dE mag ^ & 



B 3 S 



dt 



4n 



4tt 3 /V /2 



(24) 



-^mag 
Eth 

Ekin 
E v 



B x (x,z) 2 + B z (x,zf 
8vr 



dxdz, 



— p g (x,z)dxdz, 



(18) 
(19) 



-p(x, z) (v x (x, zf + v z (x, z) 2 ) da(<2@) 
J J p{x,z)g\z\dxdz, (21) 



Etot — E ma „ + E th + E kin + E„, 



(22) 



respectively. They are integrated in whole simula- 
tion region, and corrected by the energy flux at the 
top and bottom boundaries. The gravitational and 
magnetic energies are converted to the thermal and 
kinetic ones. 

The top of inflating loops collides with the anti- 
parallel magnetic field in the Galactic halo (Figs. 2b 
and 3b), so that the current sheets are created at 
t ~ 60 (Fig. 46). 

Phase III: The Magnetic Reconnection Trig- 
gered by Tearing Instability (t > 60) 

The drift velocity increases steeply at t ~ 60 (Fig. 
5 a). The tearing instability occurs at the dissipa- 
tion region, and creates the small magnetic islands 
(Tanuma 2000; Tanuma et al. 2001a). The islands are 
ejected along the current sheet. Just after the plas- 
moid ejection (t ~ 62), the drift velocity increases 
above the threshold of anomalous resistivity. The 
anomalous resistivity sets in, and Petschek-like recon- 
nection starts. The gas is accelerated horizontally to 
~ 8 — 13 (Fig. 56), which is equal to the local Alfven 
velocity. The gas is heated to 125 (Figs. 2b 

and 2 c). It is explained by 



T max ~ 



rih 



1 

PJ n out 



T- 

± in> 



(23) 



where v- m (= ev\), e, and S are the inflow velocity, 
reconnection rate, and reconnection region size, re- 
spectively. 

3.2. The Parameter Survey 

3.2.1. The Dependence on the Existence and 
Direction of the Magnetic Field in the 

Galactic Halo 

We examine the effect of the magnetic field in 
the Galactic halo, by comparing the anti-parallel- 
field model (typical model; model Al) with the no- 
magnetic-field model (model B), parallel-magnetic- 
field model (model C), and model D (anti-parallel 
magnetic field in a Galactic halo and parallel mag- 
netic field in the other Galactic halo) (Table 1). In 
models B and C, no magnetic reconnection occurs. 
In model C, the parallel magnetic field in the Galac- 
tic halo suppresses Parker instability. The velocity is 
~ 3 — 6, which is larger than that of model B (Fig. 
5), because the gas is compressed between the mag- 
netic loop and ambient magnetic field. In model D, 
the magnetic energy release rate is very small. On the 
other hand, the distribution of density do not depend 
on the magnetic reconnection. 

3.2.2. The Dependence on the Magnetic Field 

Strength 

We examine the dependence of results on B x q[= 
(87rp 9 o//?o) 1/2 ], i-e., 0o (models A2 and Fl-9; Table 
1). Figure 9a shows the /^-dependence of magnetic 
energy release rate. It is determined by the Poynting 
flux entering the reconnection region, 



dE- 



mag 



dt 



d2 

2— Sv i: 

4-7T 



cx /3 



-3/2 



(25) 
(26) 
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where S is the reconnection region size, V{ n [= ev a = 
eB/(4np) 1 ^ 2 ] is the inflow velocity to the reconnec- 
tion region, and e is the reconnection rate. It is also 

shown (~ 4.5/3 ), by assuming e ~ 0.1, S ~ 20, 
and p g ~ 0.6. It explains the results well. Figure 
96 shows the time when maximum heating rate is at- 
tained. It is determined by the time scale of Parker 

/ 1 1 / 2 

instability, i.e., Alfven time (ta oc oc /3 ' ), as 

shown (~ 100w^ ~ 90/3 ), where we assume that 
the vertical of loop is ~ O.Iwa (Matsumoto et al. 
1998) and the reconnection occurs at z ~ 10. 

Figure 9c shows the maximum temperature (T max ). 
The gas is heated to 



1 + 



1 



PoJ n ou t 



oc 1 + 



(27) 
(28) 



(see eq. [23]), as also plotted (~ 25[1 + Pq 1 ]), by as- 
suming n- Ui ~ n out and T Ui ~ T c ~ 25. It explains the 
results well. Figure 9c? shows the maximum velocity. 
The gas is accelerated to Alfven velocity, 



B 



(47rp)V2 



oc (3, 



-1/2 



(29) 
(30) 



as shown by assuming T ~ 25. It also explains the 
results well. 

3.2.3. The Dependence on the Magnetic Field in the 

Galactic Halo 

We examine the position (z^) of anti-parallel mag- 
netic field (models A2 and Gl-3; Table 1). In the 
near- magnetic- field model (model G3), the tearing in- 
stability occurs between the Galactic disk and halo. 
The magnetic field dissipates in an early phase {t ~ 
— 50), and the wavelength of the magnetic loops is 
shorter than that of the other models. 

We also examine the existence of magnetic field in 
the Galactic halo, for the near-magnetic-field models 
(models G3, K, and L; Table 1). In models K (the no- 
magnetic-field model) and G3, the magnetic energy is 
released gradually by the tearing instability or mag- 
netic dissipation in an early phase (t ~ — 35) (Figs. 
10a and 106). In model L (the parallel-magnetic-field 
model), the dissipation is suppressed (Fig. 10c). 

We, furthermore, examine (3$ in the near-magnetic- 
field models (models G3 and J 1-8; Table 1), as also 
shown in Figure 9 by A. The /^-dependence of these 
models are equal to that of high-magnetic-field mod- 
els. 



3.2.4. The Dependence on the Resistivity Model 

We examine the resistivity model. Figure 11 shows 
the reconnection rate, in the anomalous (model Al) 
and the uniform resistivity model (models M and N). 
In the uniform resistivity model, Sweet-Parker re- 
connection occurs, so that the reconnection rate is 
much smaller than that in the anomalous resistivity 
model. The reconnection rate increases transiently 
at t ~ 40. It is different from Yokoyama & Shibata 
(1996) 's results: their reconnection rate is not time- 
dependent. The reason is as follows; In our model, 
the position of magnetic field in the halo is lower than 
theirs, so that the magnetic field is stronger and the 
current density is larger at current sheet. The recon- 
nection rate, however, stays at nearly steady small 
value, which is consistent with the Sweet-Parker scal- 
ing: rj\J\ ocRe™^ 2 oc r\ l l 2 . 

Table 3 shows the route to fast reconnection. In 
the anomalous resistivity model, Petschek type re- 
connection occurs after tearing instability and onsets 
of anomalous resistivity, while Sweet-Parker type re- 
connection occurs in the uniform resistivity model. 
Figure 12 displays the dependence on the resistivity 
models. It shows a characteristic patterns of two re- 
connection models. 

4. DISCUSSION 

4.1. The Origin of X-Ray Gas in the Galactic Halo 

Parker instability is initiated by a small pertur- 
bation, a supernova explosion (Kamaya et al. 1997; 
Steinacker & Shchekinov 2001), collision of the high- 
velocity clouds (HVCs; Santillan et al. 2000), and 
cosmic rays etc. The X-ray gas can be generated, if 
the magnetic reconnection is triggered by Parker in- 
stability or the collision of HVCs (Kerp et al. 1994, 
1996). 

The reconnection heats the gas to 



T - 10 6 



n 



B 



K. 



,10- 3 cm- 3 / \3/^G, 
The reconnection also accelerates the gas to 



(31) 



v A ~ 300 



n 



1/2/ B 



lO- 3 cm- 3 / V3 liG 
The duration of fast reconnection is 
I 



t 



ev A 



3000 



kms" 1 . (32) 



(33) 



I 



-l 



100 pcj vo.i, 

X 1 / 2 / B 



10 3 cm -3 



3 fiG 



yr, (34) 
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where the reconnection rate is e(= v- m /vA). 

The anomalous resistivity sets in at least when the 
current sheet thickness becomes comparable with ion 
Lamor radius [~3x 10 7 (T/10 4 K) 1 / 2 (5/3/i G)" 1 
cm]. The radius of magnetic island (i.e., plasmoid) 
is larger than it. The island in 2D simulation is heli- 
cally twisted magnetic tube in 3D. Its volume, mass, 
and luminosity are 



V t 



tube 



10 b 



10 pcy vioopc 



P \ 3 
1 cm , 



(35) 



M tube ~ 2 x 10 d 



33 / 'P 



10 pc/ \100 pc/ \1 cm 3 



n 



L tube ~ n 2 A(T)V tube 



10 



31 



10 3 cm 3 



A(T) 



10 21 erg cm 3 s 



3 -i 



10 pc/ VlOOpc 



(37) 



where r p and Z p are radius and length of the plasmoid, 
and A(T) is the cooling function. Then, the number 
of magnetic tubes at such size is ~ 10 7 ~ 8 . 
Total released magnetic energy is 



£mag ~ 10 



51 



B \ 2 ( A 



3nGJ \lkpc) U00 



I 



pc 



erg, 



(38) 

where A 2 ot is the total area of many reconnection re- 
gions, and I are the typical thickness of magnetic loop. 
The energy release rate, then, is 



dR 



mag ^ 1Q 40 



B \ 3 
3 nGj 



\ N 2 

'Hot 



n 



1 kpc / \ 10 3 cm 3 



-1/2 



erg s09) 



It is also derived from eq. (26). It can explain 
the X-ray luminosity (~ io 39 ~ 40 erg s -1 ; Pietz et 
al. 1998). The heated gas is confined for r con( j ~ 
10 9 (n/10~ 3 cm- 3 )(A cfr /3 kpc) 2 (T/10 6 K)~ 5 / 2 yr (eq. 
[17]), where \ c g is the effective length of the helical 
magnetic tubes. The reconnection creates the bipo- 
lar jet, forming the high velocity gas at Alfven veloc- 
ity (~ 300 km s" 1 ; see also Nitta et al. 2001, Nitta, 
Tanuma, & Maezawa 2002). The Doppler shift of the 
bipolar jet will be the evidence of the Galactic flare, 
proposed in this paper. 

The plasmoid of cool gas is also created at the 
same time by the reconnection (Yokoyama & Shibata 
1996). It is also confined by magnetic field. The high 
velocity cool gas as well as hot gas will be the evi- 
dence of our model. 



4.2. The Comparison With Other Numerical 
Simulations 

The 2D model examined in this paper is an exten- 
sion from the supernova-shock driven reconnection 
model in the Galactic disk (Tanuma et al. 1999a, 
1999b, 2001a; Tanuma 2000). Some 2D numerical 
simulations were done for the reconnection in the 
Galaxy (Zimmer et al. 1997; Birk et al. 1998; Tanuma 
et al. 1999a, 2001a). They assume the current sheet 
for the initial condition (see also Ugai 1992; Ugai & 
Kondoh 2001). Different from their models, Parker 
instability creates the current sheets spontaneously 
in our present model. 

In our results, the tearing instability triggers 
Petschek reconnection. Sweet-Parker current sheet, 
however, would be created, and the secondary tear- 
ing instability will occur in the current sheet be- 
fore Petschek reconnection, if we use fine grids like 
Tanuma et al. (1999a, 2001a)s' model (Table 3). Fur- 
thermore, in the actual Galaxy (> 10 15 ), as well as in 
solar corona (~ 10 12 ), the magnetic Reynolds num- 
ber is much larger than that of the numerical sim- 
ulations (~ 10 4 ~ 5 ), so that the current-sheet thick- 
ness must become much smaller than the current- 
sheet length, to set in the anomalous resistivity (e.g., 
through "fractal tearing instability"; Tanuma 2000; 
Shibata & Tanuma 2001; Tanuma et al. 2001a). 

A different situation from the solar flare model 
is the growth of odd-mode of Parker instability. 
Another different result from Yokoyama & Shibata 
(1996)'s model is the time variation of reconnec- 
tion rate in the uniform resistivity model (studied 
in section 3.2.4). On the other hand, in our model, 
Petschek reconnection occurs in the anomalous re- 
sistivity model, while Sweet-Parker reconnection oc- 
curs in the uniform resistivity model. This result is 
the same with Yokoyama & Shibata (1994). We con- 
firmed this in more general situation. 

4.3. Turbulence and 3D Effects 

The MHD turbulence may be also important in re- 
connection problem (Lazarian & Vishniac 1999). The 
effective reconnection rate increases if the diffusion 
region is in a state of MHD turbulence. It is, how- 
ever, difficult to resolve MHD turbulence (even if in 
2D). 

Recently, we revealed that the fast reconnection oc- 
curs after the current sheet thinning by the tearing 
instability in both 2D (Tanuma et al. 1999b) and 3D 
models (Tanuma 2000; Tanuma et al. 2001b) under 
a low spatial resolution and the assumption of initial 
current sheet. We found no difference between 2D 
and 3D models in reconnection rate, inflow velocity, 
velocity of reconnection jet, temperature of heated 
gas, slow shock formation accompanied with Petschek 
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reconnection, fast shock formation due to the colli- tially supported by Japan Science and Technology 
sion between the reconnection jet and high pressure Cooperation (ACT-JST). 
gas, and time scale of these phenomena. These re- 
sults do not have a large quantitative difference from 
2D models with a high spatial resolution examined 
by Tanuma et al. (2001a) (see also Ugai & Shimizu 
1996). Rayleigh- Taylor instability is, however, ex- 
cited due to the collision, which occurs much after the 
onset of anomalous resistivity (Tanuma et al. 2002; 
in prep, for study of details). 

We can not examine 3D model with the same small 
grid with that of 2D model, so that we examine 2D 
model with a small grid in this paper. The secondary 
tearing instability, however, may be different between 
2D and 3D, if we assume enough small grid. It is also 
important to study details of 3D structure of diffu- 
sion region, the tearing instability, plasmoid ejection, 
heating of X-ray gas in the Galactic halo, and cre- 
ation of high velocity "bipolar jet", by performing 
3D simulations with enough fine grid. It is our future 
work. 

5. SUMMARY 

We propose the Galactic flare model for the origin 
of the X-ray gas in the Galactic halo. We examine the 
magnetic reconnection triggered by Parker instability 
in the Galaxy, by performing the two-dimensional re- 
sistive MHD numerical simulations. 

At the initial condition, we assume the horizontal 
flux sheet in the Galactic disk, and the anti-parallel 
magnetic field in the Galactic halo. The magnetic 
field inflates toward the Galactic halo, by Parker in- 
stability. It collides with the anti-parallel magnetic 
field in the Galactic halo, and creates the current 
sheets. The tearing instability occurs in the current 
sheet, and creates magnetic islands. Just after the 
plasmoid ejection, the anomalous resistivity sets in, 
and Petschek reconnection occurs. In the Galactic 
halo, the magnetic field of B ~ several (iG can heat 
the gas to T ~ 10 6 K. If the bipolar jet of high- 
velocity hot gas or cool gas at Alfven velocity (~ 300 
km s _1 ) is observed, it is the evidence of the magnetic 
reconnection model in the Galactic halo. 

The authors thank R. Matsumoto in Chiba Univer- 
sity and K. Makishima in Tokyo University for fruit- 
ful discussions. The authors also gratefully acknowl- 
edge the constructive advice and useful comments 
of our referee toward improving our manuscript. 
The numerical computations were carried out on 
VPP5000 at the Astronomical Data Analysis Cen- 
ter of the National Astronomical Observatory, Japan, 
which is an inter-university research institute of as- 
tronomy operated by Ministry of Education, Culture, 
Sports, Science and Technology. This work was par- 
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Table 1 

Parameters. The typical model is named Al. The unit of length is ~ 100 pc. 



model 


/3o a 








Bi c 


Ax f x 


Az s 


N x h x 




Lj 


x 




n l 


Al 


0.200 


5.0 


6.0 


n 


TT 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


anomalous 


A2 


0.200 


5.0 


6.0 


n 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


B 


0.200 


— 


6.0 


none 


none 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


anomalous 


C 


0.200 


5.0 


6.0 


TT 


TT 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


anomalous 


D 


0.200 


5.0 


6.0 


TT 


TT 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


anomalous 


Fl 


0.100 


5.0 


6.0 


1 4- 


1 4- 


350 x 


0.200 


203 x 


303 


70 




60.0 


CI VI CITY! CI 1Q 


F2 


0.175 


5.0 


6.0 


u 


u 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


F3 


0.225 


5.0 


6.0 


u 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


F4 


0.250 


5.0 


6.0 


u 


TT 


0.350 x 


0.200 


203 x 


303 


70 


X 


60.0 


anomalous 


F5 


0.275 


5.0 


6.0 


n 


U 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


F6 


0.300 


5.0 


6.0 


n 


U 


0.350 x 


0.200 


203 x 


303 


70 


X 


60.0 


anomalous 


F7 


0.500 


5.0 


6.0 


n 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


F8 


1.000 


5.0 


6.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


F9 


1.200 


5.0 


6.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


Gl 


0.200 


6.0 


7.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


G2 


0.200 


4.5 


5.5 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


G3 


0.200 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


Jl 


0.175 


4.0 


5.0 


ii 


U 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J2 


0.180 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J3 


0.225 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J4 


0.250 


4.0 


5.0 


ii 


U 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J5 


0.275 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J6 


0.300 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J7 


0.500 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


J8 


1.000 


4.0 


5.0 


11 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


K 


0.200 




5.0 


none 


none 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


L 


0.200 


4.0 


5.0 


TT 


TT 


0.350 x 


0.200 


203 x 


303 


70 


x 


60.0 


anomalous 


M 


0.200 


5.0 


5.0 


TT 


TT 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


uniform(?j = 0.05) 


N 


0.200 


5.0 


5.0 


TT 


TT 


0.300 x 


0.075 


403 x 


603 


120 


x 


69.2 


uniform(r/ = 0.10) 


a The ratio of the gas to magnetic pressure 



b The position of magnetic field in the Galactic halo 

c Thc position of hot gas in the Galactic halo [z c = Zh + 1) 

d The magnetic field in the upper corona (TT :parallel, ||:anti-parallel) 

c The magnetic field in the lower corona (|t:parallel, f|:anti-parallel) 

f The grid size in x-axis (uniform) 

g Minimum grid size in z-axis 

h Thc number of grids in x-axis 

'The number of grids in z-axis 

JThe simulation region size in x-axis 

k The simulation region size in z-axis 

'The resistivity model 
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Table 2 

Variables at the initial condition for the typical model (model Al) 



variable 


cool dense disk 






hot rarefied halo 






T 


T Q = 1 


(\z\ 


< Zc) 


T = 25T 


(\z\ 


> Zc) 


P 


Po = 1 


(\z 


= 0) 


p ~ 0.004p 


(\z\ 


~ 2ft) 


Pg 


PgO = 0.6 


(\z 


= 0) 


Pg ~ 0.lp g0 


(\z\ 


~ Zft) 


B x 


B x0 ~ 8.68 


(\z 


= 0) 


B x ~ 0.33,0 


(\z\ 


~ Zfc) 


Plot 


Ptot.o = 3.6 


{\z 


= 0) 


Ptot ~ O.lptot.O 


(\z\ 


~ Zh) 


A) 


00 = 0.2 




= 0) 


A) = 0.2 




>z h ) 
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Table 3 

The route to the fast reconnection. Parker instability trigger model is examined in this 
paper. In this case, resistivity model determines the route. The supernova shock trigger 

MODEL IS EXAMINED IN TANUMA ET AL. (2001 a). 
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Fig. 1. — Schematic illustration of our numerical simulations, (a) As the initial condition, we assume the horizontal, parallel 
magnetic field and cool, dense gas in the Galactic disk, and the nearly uniform magnetic field which is anti-parallel to the disk 
field and hot, rarefied gas in the Galactic halo. The anti-parallel field is created, for example, by Coliori's force, (b) Parker 
instability is initiated by the random perturbations in the Galactic disk. The magnetic field in the Galactic disk bends across the 
equatorial plane, (c) The inflating magnetic loop collides with the anti-parallel magnetic field in the Galactic halo. The magnetic 
reconnection occurs, (d) The heated gas is confined by the magnetic field in the Galactic halo for a long time. 

Fig. 2. — The temperature (T) in the unit of To ~ fO 4 K, magnetic field lines, and velocity vectors, for a typical model 
(model Af), The units of length, velocity, and time are f kpc, fO km s _1 , fO 7 yr, respectively, (a) Parker instability occurs in the 
Galactic disk. The magnetic field of the Galactic disk inflates to the Galactic halo, (b) The inflating magnetic loop collides with 
anti-parallel magnetic field in the Galactic halo. The current sheets are created, (c) The magnetic reconnection occurs at t ~ 62, 
and heats the gas to T max ~ f25. (d) The heated gas is confined by the magnetic field in the Galactic halo. 

Fig. 3. — Same with Figure 2, but the density (p) in the unit of I0~ 25 g. (b) The gas, which is initially between the Galactic 
disk and halo, are raised by the inflating magnetic loops, (d) The large high-density clouds are created at the valleys of magnetic 
loops, and will become star forming regions. 

Fig. 4. — Same with Figure 2, but the current density (J), (b) Eight magnetic loops collide with anti-parallel magnetic field in 
the Galactic halo, so that the current density increases there, (d) The magnetic field confines the high-current-density regions at 
the valleys of magnetic loops near the Galactic plane. 

Fig. 5. — Time variations of (a) drift velocity (v d = \J\/p) and (b) velocity (|wa:|), for model Af (solid lines; typical model, 
i.e., anti-parallel-magnetic-field model), model B (dotted lines; the no-magnetic-field model), model C (dashed lines; the parallel- 
magnetic-field model), and model D (dashed and dotted lines; the anti-parallel magnetic field in a halo and parallel magnetic field 
in the other halo), (a) In models Af and D, the drift velocity increases steeply above the threshold (v c = 400) of anomalous 
resistivity at t ~ 62 — 65. (b) In all models, the velocity increases gradually by Parker instability to ~ 2.5 — 3 (t ~ 40 — 62). In 
models Af and D, by the magnetic reconnection (t > 62), the velocity increases steeply to ~ 8 — 13, which is equal to the local 
Alfven velocity in the Galactic halo. 

Fig. 6. — Time variations of the various energies, (a) The gravitational energy is released by Parker instability in its nonlinear 
phase (t ~ 40 — 62). The magnetic energy is released mainly by the magnetic reconnection (t > 62). (b) The magnetic energy 
release rate increases to —dE ma . s /dt ~ 75 by the magnetic reconnection (t ~ 70). 

Fig. 7. — The profile of the variables in x — —13.8 at t = 70. The magnetic reconnection occurs around (x,z) ~ (—26,7.8). 
The fast reconnection such as Petschek model is accompanied by two slow shock regions along the reconnection jets. The profile 
of current density (|J|) has two peaks at two slow shock regions. 

Fig. 8. — Same with Figure7, but z = 6.5. The gas is accelerated, by magnetic tension force, to the local Alfven velocity in the 
Galactic halo, and collides with the ambient gas. The fast shock region is created around x ~ —10. 



Fig. 9. — The dependence of the results on /3o, i.e., the magnetic field strength [B x = (8ttp 9 //3o) 1/2 ] (models A2, Ff-0, G3, 
and Jf-8). The results of the far- magnetic-field models (models A2 and Ff-9) and the near- magnetic-field models (models G3 and 
JI-8) are shown by o (with the dashed lines) and A (with the solid lines), respectively, (a) The maximum magnetic energy release 
rate (— dEmag/cft)- It is determined by Poynting flux (oc /3q 3 ^ 2 ; the dashed line), (b) The time when the maximum magnetic 
energy release rate is attained. The time scales of these phenomena depend on Alfven time (oc v^ 1 oc f3^ 2 ; the dashed lines), (c) 
The maximum temperature (T max ) of heated gas. It is determined by the released magnetic energy (oc I + 1/Po', the dashed line), 
(d) Maximum v x , which is determined by Alfven velocity (va oc Pq 1 ^ 2 ; the dashed line). 

Fig. 10. — Time variations of the various energies, for the near-magnetic-field models (models G3, K, and L). (a) In model G3 
(the anti-parallel-magnetic- field model), the magnetic field dissipates by the tearing instability or magnetic dissipation (t ~ — 40). 
The wavelength of the magnetic loops is shorter than that of the typical model, (b) In model K (the no-magnetic-field model), 
the magnetic dissipation occurs in an early phase, (c) In model L (the parallel-magnetic-field model), the magnetic field does not 
dissipate in an early phase, because the parallel field of the Galactic halo suppresses the dissipation. 

Fig. If. — The time variation of the reconnection rate (»y|J|), in the anomalous resistivity model (model Af) and the uniform 
resistivity model (rj = 0.05 [the dashed line; model M] and r\ — 0.10 [the dotted line; model N]). Petschek reconnection occurs in 
the anomalous resistivity model. On the other hand, Sweet-Parker reconnection occurs, in the uniform resistivity model. 

Fig. 12. — The current density, magnetic field lines, and velocity vectors, for (a) the anomalous resistivity model (model Af) 
and (b) the uniform resistivity model (r) = 0.05; model M]). (a) Petschek reconnection occurs in the anomalous resistivity model. 
The small dissipation region and slow shock regions are the characteristics of Petschek reconnection. (b) In the uniform resistivity 
model, the long current sheet forms because Sweet-Parker reconnection occurs. 
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